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Abstract 

Determinations of the size and shape of the interaction region from 
fc-particle (k = 1,2,...) momentum distributions of identical particles 
are analyzed. The full group of transformations changing the single 
particle density matrix without affecting any of the measurable mo- 
mentum distributions is identified. The corresponding uncertainties 
in the inferred parameters of the interaction region are discussed. 

PACS numbers 25.75. Gz, 13.65.+i 
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1 Introduction 

The study of Bose-Einstein correlations, or more generally of single particle, 
two particle etc. momentum distributions for identical particles, is the stan- 
dard tool used to learn about the size and shape of the interaction region 
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where the hadrons are produced. We will explicitly discuss spin zero bosons, 
because in practice this is the most important case. Much of the argument, 
however, can be extended to bosons with non-zero spins and to fermions. 

Most published analyses are variants of the following approach which we 
will also follow. One defines a single particle density matrix in the momen- 
tum representation p(p;p'). The diagonal elements of the (symmetrized) 
72,-particle density matrix are taken as PQ 

n 

p(Pl,---,Pn;Pl>--- 5 Pn) =^2Ylp(PpPPj)^ W 

P j=l 

where the summation is over all the permutations j — *> Pj of the indices j. 
The n-particle momentum distribution is given by the diagonal elements of 
this ri-particle density matrix 1 . Thus, all the information about the system 
is contained in the single-particle density matrix p(p;p'). The question is: 
how well this function can be determined from the experimental data? 
The single particle momentum distribution is 

P( P ) = p(p;p). (2) 
Thus the diagonal elements are measurable. The two-particle distribution is 

p (pi, p 2 ) = p(pi; PiMp 2 ; P2) + Ip(pi; p 2 )I 2 , (3) 

where the hermiticity of the density matrix 

p(pi; p 2 ) = p*(p 2 ; Pi) (4) 

has been used. Thus, the absolute value of every matrix element is also 
measurable. For the three-particle momentu distribution formula yields 

P(Pi; P 2 , p 3 ) = . . . + 23H [p(Pi; p 2 )p(p 2 ; P 3 )p(P 3 ; Pi)] , ( 5 ) 

where the dots denote terms which can be determined from the single par- 
ticle and two-particle momentum distributions and 3? stands for real part 

^^One can do somewhat better by including more-particle density matrices, but these 
corrections, known as multiparticle corrections, are usually neglected, because they are 
not very important (cf. e.g. the review [2]). 
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of. In general, since every permutation can be decomposed into cycles, the 
n-particle momentum distribution is 

^(Pi, • • • . Pn) = • • • + 23? [p(Pi, p 2 )p(p 2 ; Pa) • • • P(Pn! Pi)] , (6) 

where the dots denote terms which can be determined by measuring momen- 
tum distributions for less than n particle. 

In a previous paper [2] it has been pointed out that the measurable quan- 
tities do not change when the substitution 

P(P5 P') - e^-f^p{ V] p') (7) 

is made 2 . In the present paper, we show that this group of transformations 
supplemented by the transformation 

p(p;p')^p*(p;p') = p(p';p) (8) 

generates the full group of transformations which leave the measurable dis- 
tributions invariant. 

We will discuss the following problem. Somebody found that a density 
matrix Po(p;p') gives the best fit to all the measured momentum distribu- 
tions. What other density matrices p(p; p') give exactly the same fit and how 
the conclusions concerning the interaction region are affected when po(p; p') 
is replaced by p(p;p')? This is the problem of the theoretical uncertainty 
in the determination of the parameters of the interaction region from the 
measured momentum distributions. 



2 What can and what cannot be measured 

Let us write the single particle density matrix in the form 

p(Pi;p 2 ) = |p(Pi;p 2 )|e ix(p - P2) . (9) 

The modulus is measurable, thus all the ambiguities result from the phase. 
It is convenient to choose for the phase the convention 

-7r<x(p,p') <tt. (10) 

2 This remains valid also when the multiparticle corrections are included. 
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The hermiticity of the single particle density matrix implies that 

x(p,p') = -x(p',p)- (11) 

Moreover, the phases x(PjtPk) are constrained by the measurable three- 
particle momentum distribution which yields 

cos [x(pi; p 2 ) + x(p 2 ; Ps) + x(p 3 > Pi)] • (i 2 ) 

The sign of the argument of the cosine is not constrained. This implies that 
p(p; p') can be replaced by its complex conjugate without affecting the fit to 
the data. 

Once the sign is fixed, it is possible to determine from the data the func- 
tion 

*(pi, P2, p 3 ) = x(pi; P2) + x(P2; p 3 ) + x(p 3 , Pi)- (is) 

Formally, this is an inhomogeneous, linear functional equation for the func- 
tion x(P)P')- us assume that some function Xo(p>p') is a solution of 
this equation. If the model agrees with experiment such a function must 
exist. Let us define further a function Xh(PtP') which is a solution of the 
homogenous equation 

Xa(PiI Pa) + Xh(p 2 ; Ps) + Xh(p 3 , Pi) = 0. (14) 

Function Xo(p ; p') + Xh(p, p') is then also a solution of equation (JT3J). Thus, 
the question what is the uncertainty of solution xo reduces to the question 
what is the full set of solutions of equation (JHJ). 

Differentiating both sides of this equation with respect to p±j and p 2 k we 
get the equations 

%%^ = 0, j,k = «,»,,. (15) 
The general solution of this equation system can be written in the form 

X(P,P') = /(P)-S(P'), (16) 

where /(p) and g(p) are arbitrary functions. Every solution of equations 
(JHJ) is a solution of equation (|T5|). but not the other way round. Substituting 
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solution ((inj) into equation (|14p one easily finds that the equation is satisfied 
if and only if 

g(p) = f(p)- (17) 

The same conclusion follows from the symmetry Thus, the most gen- 

eral single particle density matrix giving the same fit to the data as matrix 
Po(P) p') can be written in the form 

p(p;p') = e i (/(Pi)-/(P2)) /0o (p;p') or e- i(/(Pl) - /(P2)) Po(p;p), (18) 

where /(p) is an arbitrary real function of the momentum p. The four- and 
more particle distributions do not introduce any more constraints. 

To see how this ambiguity can be eliminated by a theoretical assumption 
let us suppose that there are no position-momentum correlations. Then the 
Wigner function (cf. (J36j0 is a product of a function of x only and a function 
of K only. Consequently, each of the density matrices p(p; p') and po(p; p') 
is a product of a function of p + p' only and a function of p — p' only. This 
implies that 

/(p)-/(p') = a-(p-p'), (19) 

where a is an arbitrary constant vector. The ambiguities in the interaction 
region, corresponding to the transformations (|18|). reduce to rigid translations 
by the vector a and/or space inversion. Since these do not affect the size and 
shape of the interaction region, one may claim that there is no significant 
ambiguity left. 

We conclude that, if somebody fitting the experimental data and using 
some model (this is necessary to get a well-defined result), finds the effective 
single particle density matrix po(p; p')> another model leading to any of the 
density matrices (fTHj) agrees with the data just as well. Of course, a priori 
some model may be more plausible than another, but these are arguments 
on a different level - unrelated to the fitting of momentum distributions. 

The effects of this uncertainty on the interpretation of the data can be 
spectacular. It is plausible, for instance, that the momenta and positions at 
birth of the hadrons are correlated so that [I], [3] 

p « Ax, (20) 
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where A is a positive constant. Using the relation between the density matrix 
and the Wigner function one easily finds that the complex conjugation of the 
density matrix corresponds to the change of sign of the coordinate in the 
Wigner function. Using further the interpretation of the Wigner function 
as an approximation to the phase space distribution one finds that models 
where 

p « -Ax (21) 
fit the data exactly as well as the previous ones. 

3 Moments of the particle density distribu- 
tion in the interaction region 

The density of particles in the interaction region is given by the diagonal 
elements of the effective single particle density matrix in the coordinate rep- 
resentation: 

p(x;x) = p(x) = /^f^>(K,q), (22) 

where 

K=i(p + p'), q = P"P' (23) 

and p(K, q) stands for p(p;p'). Integrating n times by parts over the com- 
ponents of q and assuming that p(K, q) and all its derivatives with respect 
to components of q tend to zero for |q| — > oo one finds in particular 

« x ^§/w e ' qx ($ p(K ' q) )' '' = *•»•'■ (24) 

Multiplying both sides of this identity by x" and integrating over x one gets 
the moments 

J dx p(x)^ = (x«> = i n jdK (|^P(K, q)) . (25) 
Usually only the first and second order moments are considered. 
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Let us introduce the following notation for the moments corresponding 
to the density matrix po(p; p')- Putting 



i JdK \J^-po(K, q)^ , j = x,y,z; (26) 



q=0 

d 2 



= -J dK {d^r k po{K ^)^ (27) 

and assuming that the effective single particle density matrix po(Pi> i P2) i s 
the correct one, the moments are 

(x) = r (28) 
( Xj x k ) = R 2 0jk . (29) 

It has been known for a long time (cf. e.g. [2]) that it is not possible to 
determine either (x) or (xjXk) from momentum measurements. The usual 
way out (cf. e.g. [2]) is to calculate instead the variances 

o*( Xj ) = (x 2 ) - ( Xj ) 2 . (30) 

Another set of much discussed quantities is related to the integrands of 
the integrals given above 



(4-p»(K.q: 



a- 



2 



;Po(K,q; 



v dqjdgi. I - 1 

^< K > = ~ P„(K,0) < 32 » 
The corresponding variances and correlation functions are 

^( K )-r 0j (K)r ofc (K) = -(^-logp (K,^ ^. (33) 

Denoting the left-hand side of this equality by R 2 HBT j k , and noting that 
according to all the second derivatives of the phase with respect to 
components of q vanish at q = 0, we get 

flW^KH-^logMi^l)^. (34) 
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It is seen that these functions are measurable and do not depend on the 
specific choice of the function po(K, x) among the functions giving the same 
fits. They are just the well-known HBT radii. We have kept the traditional 
notation with the squares (cf. e.g. 0) in spite of the well-known fact that 
Rh BT for j 7^ k may be negative (cf. e.g. pj). An important point is 
that, they are not sufficient to calculate the overall variances and correlation 
functions because in general 

<r 0i (K)r ofc (K)> ^ (r 0i (K))<r 0fc (K)>. (35) 

The equations for the moments become more suggestive when written in 
terms of the Wigner function, related to the density matrix po(K, q) by 



p (K,q) = Jdx WVK,x)e- tqx . (36) 

Then 

Jdx Xj W (K,x) 
r ° j(K) " fdxW (K,x) ' (37) 

^ fc(K) = /dxWb(K,x) • (38) 

Interpreting the Wigner function as the phase space density, as often done, 
one sees that rgj(K) and i^ - fc (K) are the moments of the components of x 
evaluated at fixed K. This point of views is supported by the fact that when 
averaged over K, i.e. multiplied by the distribution of K (p (K) = po(K, 0)) 
and integrated over K, they reproduce correctly the overall moments. How- 
ever, according to Heisenberg's uncertainty principle, when the variance of 
K tends to zero the variance of x tends to infinity. Only when the variance 
of K grows from zero to finite values, destructive interference between the 
contributions corresponding to different K values limits the x-range. Thus, 
the interpretation of r j(K) and i?Q - fe (K) as moments of x at given K has 
heuristic value, but should be used with care. Note also that the average 
over K of the variance cr 2 K (xj) does not give a 2 (xj). 

In the following section we will discuss how the moments introduced here 
are affected by the uncertainty in the phase of the density matrix p(p, p'). 
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4 Uncertainties in the determination of the 



moments 

Let us see how the first and second order moments of the coordinates change 
under transformations (|18|) . The terms higher order than second in the 
components of q do not contribute and, therefore, the factors multiplying 
Po(K,q) and pg(K,q) can be written respectively as 



w (tc \ i + -V^ d/( K ) 1^ g/gQ df(K) 
F ± (K,q) = 1 ±,J>_ - - £ ^ (39) 

j J j,k J 

Using the formulae from the preceding section one finds 



4(K) = ^(K)-^(K)^^ + r ofe (K)^^) + 
df(K) df(K) 



dKj dK k 



(40) 
(41) 



The HBT radii remain, of course, unchanged. 
The formulae for the overall moments are 



(42) 

/ a/(K)a/(K) \ 

\ ax, 3/4 /„' (4d) 

where (. . .} x denotes averaging over K. The overall variances and correlation 
functions are 
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-r-n - 1# (K)) + /?m?m\ ./?m\ /?m 



r (K) df(K) \ + r l 9 ^ 

9/(K)\ /df(K) 
r ok (K)—— ) +r ok ' 



K 



d Kj / K ^'° k \ dK 3 /A 

They are not measurable without further assumptions. 

Let us discuss some illustrative examples. For the models where 



ro(K) = 0, (45) 



one gets 

Si ~~ 1 l HBT,jj 



R% = C,,, + «f(Kj), J = V, z, (46) 



where the variance 




/ df(K) \ / df(K) 
\ dK 3 J K \ dK, 



(47) 



K 



It is seen that the HBT radii are a rigorous lower bounds for the true radii, 
while there is no upper bound. The interaction region can be assumed to 
be the size of a football and agreement with experiment remains as good as 
before. In practice, of course, so large interactions regions are excluded by 
what is known about the detector, but this is information going beyond the 
measurements of momentum distributions. 

If a model giving exactly the same fit as the models satisfying ()45j) exists, 
it must be possible to find a function /(K) satisfying 

V/(K) = r (K). (48) 

Unless there are singularities, this equation can be solved when the vector 
field r (K) is rotationless: 

V x r (K) = 0. (49) 



10 



Consider the following model discussed by Akkelin and Sinyukov jHj. A 
classical gas of noninteractring particles has at time t = a Gaussian dis- 
tribution in coordinate space with (x) = and a Boltzmann distribution in 
momentum space; there are no position-momentum correlations. It is obvi- 
ous how this gas will evolve. The momentum distribution will not change. 
In coordinate space each group of particle having identical momenta K will 
propagate rigidly with velocity — , where m is the particle mass. This picture 
can be confirmed by solving the corresponding Boltzmann equation j^j. For 
the particles with momentum K 

(x)(K) = -t = r (K), (50) 
m 

where the second equality corresponds to the choice /(K) = 0, and corre- 
sponds to 

Rjk = R 2 HBT,jk- (51) 

We can choose, however, just as well 



K 

f(K) = — t. (52) 



The corresponding density matrix satisfies and consequently yield; 



s 



2 



Rjj - R HBT,n + (jrf* 2 ) ■ ( 53 ) 

From the description of the model, which contains more information than 
just about the momentum distributions, we know that the first description 
is wrong and the second is correct. In this example the HBT radii can be 
interpreted as follows. For particles with any given value of K there is some 
distribution of x. Shifting rigidly all these x distributions on top of each 
other so that their averages (x)(K) coincide one gets a reduced distribution. 
The HBT radii are the radii of this reduced distribution. 

In jj] a model was proposed, where for the transverse degrees of freedom 

(xj^Kj) = XKj (54) 
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and A does not depend on the transverse momenta. Again, using the descrip- 
tion of the model one finds 

m) = \k) (55) 

and 

R% = R'hbtjj + A 2 (K*) K . (56) 

This was used to explain why experimentally the HBT radii decrease when 
the particle mass increases. In this model the true radii R 2 ^ do not depend on 
the particle mass and the observed decrease of the HBT radii with increasing 
particle mass just reflects the well-known mass dependence of the second term 
on the right-hand side. 

5 Conclusions 

We have considered the most popular class of models, where all the infor- 
mation about the size and shape of the interaction region is contained in a 
single-particle density matrix p(p; p') . The aim is to determine this matrix 
from momentum distributions measured for sets of identical particles. We 
show that, without further assumptions, this is not possible. For each single 
particle density matrix there is an equivalence class of matrices which give ex- 
actly the same fits to all the momentum distributions. All these matrices are 
related by the transformations (fTHj) . Thus, the full group of transformations 
leaving the momentum distributions invariant has been identified. 

The HBT radii, which can be unambiguously determined from the mo- 
mentum distributions, are in general different from the true radii of the inter- 
action region. Sometimes, but not always, they are lower bounds for the true 
radii. In order to determine the true radii it is necessary to determine a func- 
tion V/(K). If only the results of momentum measurements are available 
function /(K) is completely arbitrary. Models, however, usually contain in- 
formation going beyond that obtainable from momentum measurements and 
then function V/(K) may be possible to determine. It is recommended, 
when proposing a model, to formulate clearly and justify the assumptions 
which yield /(K) = 0, because these assumptions cannot be justified even 
by the best fits to the measured momentum distributions. 
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We have discussed a comparatively simple version of the standard ap- 
proach. It is easy to extend our analysis to include the multiparticle correc- 
tions. The final state interactions and, more generally, interparticle corre- 
lations, however, would make the discussion much more complicated. It is 
most unlikely, nevertheless, that they would make the result simpler. Ambi- 
guities analogous to those described here are expected to occur also in these 
more complicated models. 
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